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Abstract. Tlie theoretical treatment of Rydberg states in one-electron ions is facilitated 
by the virtual absence of the nuclear-size coirection, and fundamental constants like the 
Rydberg constant may be in the reach of planned high-precision spectroscopic experiments. 
The dominant nuclear effect that shifts transition energies among Rydberg states therefore 
is due to the nuclear mass. As a consequence, spectroscopic measurements of Rydberg 
transitions can be used in order to precisely deduce nuclear masses. A possible application 
of this approach to the hydrogen and deuterium, and hydrogen-like lithium and carbon is 
explored in detail. In order to complete the analysis, numerical and analytic calculations of 
the quantum electrodynamic (QED) self-energy remainder function for states with principal 
quantum number n = 5, . . . , 8 and with angular momentum i = n — 1 and £ = n — 2 are 
described (j = £ i: ^). 
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1. Introduction 

The celebrated high-precision spectroscopic experiments of atomic hydrogen H] |2] |3] ID |5] |6] 
|7] m |9l do) typically involve transitions among quantum states with low principal quantum 
numbers, e.g., the hydrogen 1S-2S transition. Unfortunately, the rather large uncertainty of 
the root-mean-square (RMS) charge radius of the proton limits the precision to which the 
nominal accuracy of spectroscopic experiments can be used in order to infer fundamental 
constants: while the most precisely measured transition in atomic hydrogen is known to 
an accuracy of 1.8 parts in 10^^, the most accurately determined fundamental constant (the 
Rydberg) is known only to an accuracy of 6.6 parts in 10^^ (see Ref. ITTl ). The subject 
is not without intricacies: Even the correct definition of the RMS charge radius for low- 
energy processes in atomic systems is a rather fundamental problem which necessitates a 
careful study of form factors which are a priori defined for scattering processes lfT2l [131 [141 . 
Indeed, in the 2006 adjustment of the fundamental constants 1 11 1, the RMS charge radius of 
the proton and deuteron are inferred by fitting the 23 most accurately measured transitions to 
QED theory {TSl . using the method of least squares. 

For the low-lying states canonically used for high-precision spectroscopy experiments, 
the dominant nuclear effect is from the nuclear size effect on the Dirac energy. By contrast, 
the QED theory of transitions among highly excited states of one-electron ions is simplified 
because of the virtual absence of the problematic nuclear size effect, and the dominant nuclear 
proerty that influences the spectrum is the nuclear mass. This effect is well understood, and 
therefore, the theory of Rydberg transitions among states with high angular momenta can 
be formulated very precisely. In Refs. ITSIITTI . Rydberg states have been proposed in order 
to avoid these limitations and to improve the accuracy of the Rydberg constant. If we are 
interested in determining the Rydberg constant, then it is important to consider hydrogen-like 
ions with well-known nuclear masses. Conversely, it becomes possible, using precise QED 
theory, to infer the nuclear mass from spectroscopic measurements, if that mass is not known 
to sufficient accuracy by other methods. That latter aspect forms the basis of the current 
investigation and the current proposal which is being investigated in this work. In selecting 
suitable atomic Rydberg transitions for a conceivable determination of nuclear masses by this 
method, it is important to consider frequencies accessible to optical frequency combs, as they 
constitute one of the most accurate devices available for measurements. 

Consequently, we here investigate the possibility to use high-precision spectroscopy in 
order to measure nuclear masses and electron to nucleus mass ratios using transitions among 
Rydberg states of hydrogen-like ions, with a special emphasis on the example cases of ^H, 
^H, ^Li and ^''C. These masses are of great interest. The mass ratios of and are 
interesting for the frequency comparison of the transitions of hydrogen and anti-hydrogen 
which are being pursued at the moment [TSl, as well as for many other applications of 
spectroscopy in general. ^Li is important for the study of solar neutrinos and for determining 
their masses |[T9l . There is currently an interesting discrepancy between the recently obtained 
mass measurement f20| and the value in the Atomic Mass Evaluation |21 1 (AME2003), with 
the more recent value from Ref. [ 20 1 being more precise. Thus, a measurement of the ^Li mass 
with a completely independent method would seem worthwhile. The mass of ^°C meanwhile 
plays a role in super-allowed beta decays which are studied in order to measure Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements 122. 23 J . The possibility of determining nuclear 
masses by high-precision spectroscopy has previously been mentioned in Refs. [l24l [25l . in 
the context of molecular spectrosocpy where the theoretical challenges in providing accurate 
predictions for mass determination appear to be much greater 

In order to determine the mass of the nucleus, we consider essentially two options. 
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These are based on the elimination of variables in the theoretical expressions for the transition 
frequencies and can be summarized as follows. Namely, up to very high accuracy lO^^"*), 
the Rydberg transition frequencies in one-electron ions are determined by the following 
parameters: (i) the nuclear mass and (ii) the Rydberg constant. Our proposed Method I 
applies if at least one isotope of a given charge number has a very well known mass (of 
relative accuracy ^ 10^^") which can be used as a reference for the determination of the 
masses of other isotopes. We then propose to measure the same Rydberg transition in 
two ionic systems corresponding to two different nuclear isotopes, one of which acts as a 
reference, and to solve the system of the two equations for two variables, namely the Rydberg 
constant and the unknown nuclear mass. Method II applies if we suppose that recent efforts 
in measuring an improved value of the Rydberg constant (see Refs. 1161 [1711261 ) are crowned 
with success, and that an improved value of the Rydberg constant (of relative accuracy 
10"^'^ . . . 10"^"') is available. In that case, the repetition of a measurement of the same or of 
a different Rydberg transition in a one-electron ion of a different nuclear isotope or even with 
a nucleus of a different charge number can directly lead to a determination of the mass of that 
nucleus (we essentially solve two observational equations corresponding to the two measured 
frequencies for the two unknowns, namely, the Rydberg constant and the nuclear mass to be 
determined). Method III is essentially identical to method II with the only difference being 
that the theoretical expression is solved for the electron to nucleus mass ratio instead of the 
nuclear mass. 

So far, transitions in circular Rydberg states have been measured up to a relative accuracy 
of 2.1 X 10^^^ in an 80 K atomic beam of hydrogen lIZTl l28l . The measurement was 
carried out in the millimeter region and shows that it is in principle possible to conduct high- 
precision measurements on Rydberg states. There is no fundamental obstacle to improving 
the experimental accuracy of infra-red and near-optical transitions among Rydberg states 
of hydrogen-like ions up to the level of lO^^"* . . . 10^^^. Using femto-second lasers and 
frequency combs, this level of accuracy has already been attained for the lower-lying states 
of hydrogen. Spectroscopy on this level of accuracy has the potential to assist in the 
determination of nuclear masses, as we show here. In order to carry out the current study, 
a supplementary calculation of the self-energy remainder function is carried out in the range 
of low nuclear charge numbers Z = 3, 6, 8, for excited states with principal quantum number 
n = 5, . . . , 8 and orbital angular momenta I = n — 1 and I = n — 2. 

The outline of this paper implies a discussion of QED theory in Sec. |2] with a detailed 
discussion of the proposal for the nuclear mass determination following in Sec. |3] The 
numerical and analytic methods for the treatment of the self-energy remainder functions are 
described in Sec. 12.21 Numerical examples concerning the attainable accuracy for the nuclear 
mass determination are given in Sec. 13.21 Finally, conclusions are drawn in Sec.|4] 

2. QED Calculations 

2.7. Status of Theory 

In formulating the QED theory of Rydberg transition frequencies lITTl [T6l [TtI |29l [30l . we 
write the transition frequency between two Rydberg states of a one-electron atom as 

Vl^2 ^V2-Vi (1) 

for a transition between quantum states |1) and |2). The quantities Vi (i — 1, 2) are related 
to the bound-state energies as i^i = Ei/h, where h is Planck's constant. We now discuss the 
different contributions to the bound-state frequencies z/^ one-by-one, recalling a number of 
expressions from Refs. ifTTl [T6l [iTl |29l l30l which are relevant for the current study. 
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The dominant contribution arises from the Dirac energy, which we convert to a Dirac 
frequency z^p. We subtract the rest mass and correct for the nuclear mass. For a state with 
principal quantum number n, total angular momentum j, and angular momentum £, we have 
the well-known result 



2 



2[l + r{Af)Y 



(2) 



a^mec/2h is the Rydberg constant, r{Af) — 



1 + r(AA) ' 

where c is the speed of light, i?o 
TOe/mAr(A/') is the mass ratio of the electron to the nuclear mass, and a is the fine-structure 
constant. In writing mN{Af), we use the subscript N in order to denote the nuclear mass, and 
reserve the argument J\f in order to differentiate a specific nucleus under investigation, e.g., 
J\f = ^Li or A/" = ^°C. The function f{n,j) is given as 



(Za) 



{n-j-'^+,J{j + ^r-{Za)^y 



(3) 



For non-S* states, there is a further nuclear-mass dependent contribution from the so- 
called Barker-Glover (BG) term which originates from the two-body Breit Hamiltonian and 
eliminates the {n, j) degeneracy of Dirac theory 1311 , 

RooC riMfZ^a" 



l + r{N) n3 [l + r{J\fW 



(l-<5,o) 



VKK 



[1 + r{U)Y 
7 1 



1 + r{N) TTV? 

g 

3 ' ^ 3^(^+l)(2£+l) 

+ TTZa [1 + r(M)] 

1(1+1) 2 



(4) 



The relativistic -recoil (RR) correction changes the frequency of a level with t > 2 

by nmnaiMiEa 

RooC 2r{M)Z^a^{ 1 



(5) 



n? (4/2 - l)(2^ + 3) 
where In/co is the (nonrelativistic) Bethe logarithm that depends on n and I. The QED 



radiative corrections for these levels contribute 

RooC 2Z^o?- 



1 



1 + r{M) TV' I 1 + r(M) k{21 + 1) 

1 4, , , ^. 32 

--lnfcoK^) + y 



[l+r(A/')]2 TT 

3^2 -£(£+1) (2^-2)! 
^^2 (2^ + 3)! 

+ {Zaf G{Za 



{Zaf In 



1 + r{M) 
{Zaf 



(6) 
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where k = (— (^j + i) is the Dirac angular quantum number and the electron 
magnetic moment anomaly. The semi-analytic expansion of the self-energy remainder 
function G{Za), which depends on n, I and j, reads 

G{Za) ^ Ago + {Zaf [A^ \n{Za)-^ + Ago + • • •] 

+ -B60 + ..-+ (- C^60 + .-. . (7) 
TT \7r/ 

Here, we use the commonly accepted notation for the QED correction terms. The letters 
denote the loop order, i.e., A coefficients arise from one-loop diagrams, B from two-loop, 
and C from three-loop QED corrections. Vacuum polarization is negligible for Rydberg 
states with 71 > 5, ^ > 3 and Z < 8, on the level of 10~^^ of relative accuracy for the 
transitions under study. So, we can restrict the discussion to the self-energy contributions to 
the A coefficients. The first subscript of the A coefficients denotes the power of Za, while the 
second denotes the power of the logarithm ln[(ZQ;)^-^]. To achieve an appropriate accuracy 
using this expansion, the Aqq coefficients have to be determined. The higher-order one-photon 
terms except the vacuum polarization contribution to A^q ||36l , which is very small, are so far 
unknown. Calculations also exist for iJgo (see Ref. |371), however only for states with ^ < 5. 
The three-loop term Cgo as well as the higher-order two-photon contributions not listed in 
Eq. (|7]) are unknown. 

Several results for the Ago coefficient of Rydberg states in the region n = 9, . . . , 16 
of principal quantum numbers have recently been obtained (see Refs. |fT6l [iTl [38l ). The 
calculation has to be done individually for each one of the ionic states, because several 
logarithmic sums over the spectrum of virtual excitations ("relativistic Bethe logarithms") 
can only be calculated numerically. In the current investigation, a different range of principal 
quantum number n — 5, ... 8 and orbital angular momenta I = n — 1 and £ = n — 2 will 
be investigated, because the corresponding transition frequencies match the most appropriate 
regions for the application of frequency combs to the most interesting isotopes. We investigate 
the one-electron systems of lithium, carbon and oxygen. The analytic results for ^go are 
compared to a nonperturbative (in Za) evaluation of the self-energy remainder function 
GsE [Za], which corresponds to the entire contribution of higher-order terms to the one-loop 
self-energy shift for the ionic states under investigation (see Tab. [T). The function Gse(^<^) 
is equal to the sum of the higher-order A terms, including Ago, due to the electron self-energy, 
as written in Eq. dTji. 

For nuclei with nonzero nuclear spin another contribution arises from the interaction of 
the nuclear spin / with the total angular momentum j of the electron. This leads to the known 
hyperfine shtting of the energy levels which is given as [39J 

Rooc Z'^o? r r(7V) K 

t'hfs 



l + r(A/') rfi \\ + r{M)\K 



(2k+ 1)(k2 - 



4' 



Y^-^tA^ 

where g is the nuclear g factor, F — I+j the total angular momentum of the one-electron ion, 
7 = ^ — (Za)'^ and N = {n — \k\Y + 2(n — |k|)7 + k^. The term of relative order 
a has been obtained in Ref. [40] ; the generalization to arbitrary quantum numbers is found 
here. If necessary, the corrections to the hyperfine sphtting up to relative order a{ZaY can 
be calculated using a generalization of the approach of Ref. ll4T1l . 
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Having discussed all the required corrections for the energy levels of Rydberg states in 
one-electron ions for the purposes of our investigation, we can write down the total frequency 
of a specific level as the sum 

Vt = l^B + J^BG + i^RR + t'QED ■ (9) 

where i = 1, 2 in the spirit of Eq. (|T|). We here assume that the lyhis has been subtracted 
from the transition frequency, if necessary, the hyperfine-fine structure mixing terms can be 
calculated according to the approach outlined in Sec. Ill of Ref. liTSl . However, we note that 
the state with the highest projection Alp in a given manifold of total electron+nuclear angular 
momentum F does not receive any corrections due to hyperfine-structure fine-structure 
mixing. 



2.2. Calculation of The Self-Energy Remainder 

The expression (|6]l for the QED corrections to an ionic level includes recoil corrections in 
the form of the electron-to-nuclear mass ratio r{N). For the evaluation of the self-energy 
correction, it is much more practical to use the non-recoil approximation, which implies an 
infinitely heavy nucleus and therefore the limit r{Af) 0. In the non-recoil limit, the one- 
photon self-energy shift AE{n£j ) of principal quantum number n, orbital angular momentum 
£ and total angular momentum quantum number j can be written as 

AE{n£j) = 2hRooC A40 (10) 



+ {Zaf (Aei In [(Za)-^] + Gs^iZa)) | , £>2, 

where the first subscript of the A coefficients denotes the order of Za while the second denotes 
the order of the logarithm ln[(Za)^^]. We have 

lim GsE(^a) = Ago. (11) 

Za^O 

The form of the expansion in Eq. ( flOl l is valid for states with orbital angular momentum 
quantum number t > 2. The known results for A40 and Aqi for Rydberg states read as 
follows, 

^^" = -^"^°- 2.(2^+1)' ^^'^ 

_ 32 3n^-^(£+l)(2£-2)! 
^''-^ (27T3)!- ^^^^ 

In Eq. the prefactor [1 + r{M)]^^ as well as the 1 + r{M) in the logarithm account for 
the reduced-mass effects. Also, we reemphasize that these results are exclusively due to the 
self-energy effects, and that vacuum-polarization shifts vanish for states with £ > 3 up to the 
order Q;(Za)i". 

The evaluation of Ago proceeds via a scale-separation parameter (overlapping parameter) 
that is used in order to express the self-energy shift AE{n£j ) in terms of two contributions, 
one from high-energy virtual photons and one from low-energy photons whose energy is 
commensurate with the atomic binding energy. The overlapping parameter cancels when 
both contributions are added, no matter what regularization is used |38|. Methods from 
effective field theory can be used in order to simplify the calculation of the high-energy part 
drastically ll42l l43l [38l . The method of evaluation has otherwise been described in detail 
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in Refs. Il35l [161 [TTl . For Rydberg states, it becomes almost indispensable to formulate 
the logarithmic sums over the virtual excitations, which contribute to Ago, in terms of a 
discretization of the spectrum on a lattice |l44l|. In Table [T] we present results for the states 
under investigation, which are highly excited Rydberg states with principal quantum numbers 
n = 5, . . . , 8 and with the configurations £ ^ n — 1 and £ = n — 2, and j = £±^ (four states 
for a given n). 

In columns 4-15 of Table [T] we present the results of a nonperturbative (in Za) 
calculation of the self-energy remainder function GsE{Za) for one-electron ions with nuclear 
charge numbers Z ~ 3, 6, 8 (lithium, carbon and oxygen), for the same states that are also 
relevant for the analytic calculation. Note that in view of Eq. ( fTTl i. the Aqq coefficient 
constitutes the limit as Za ^ of a quantity which is otherwise dependent on the nuclear 
charge. Therefore, in Table [1] a specific value of Z has to be given along with any Gse 
(the value of a employed in the calculation is = 137.036). For Rydberg states, the 
dependence of GsEiZa) on its argument is very weak, and the deviation \GsE{Za) — Aqq\ 
numerically is less than 10"'* for all states and all nuclear charge numbers under investigation 
here. Note that the evaluation of the self-energy remainder function for Rydberg states 
is a highly nontrivial problem [see Il45l 1. In general, when treating Rydberg states both 
analytically as well as numerically, one has to be very careful not to terminate the sum 
over the angular momenta of the virtual states too early, because numerically significant 
contributions are due to virtual excitations with high angular momenta (displaced from the 
angular momentum of the reference state by one unit). 

2.3. Estimate of Theoretical Uncertainties 

In order to gauge the applicability of our method for determining nuclear masses, we have 
to investigate the inevitable theoretical uncertainties that affect the theoretical predictions 
of Rydberg transition frequencies in hydrogen-like ions. These theoretical uncertainties 
ultimately determine the limits of accuracy for the extraction of the nuclear masses. 

Two different sources of theoretical uncertainty have to be distinguished. The first of 
these comes from the input parameters necessary to evaluate the theoretical expression in 
Sec. 12.11 For these input parameters (nuclear masses and other fundamental constants), we 
use the values of the fundamental constants from COD ATA 2006 fTTl, the masses from the 
2003 Atomic Mass Evaluation [21 1 (AME2003) and the recent measurement for the ^Li mass 
||20l which are given in Table |2] Moreover, the experimental value for is used in the 
theoretical expression (|6]l; it reads 

fle == 1.159 652 180 73(28) x 10"^ (14) 

as obtained with a one-electron cyclotron [46 1. We also remember that the 2006 CODATA 
value of the Rydberg constant (see Table|2]i carries a relative accuracy of 6.6 x 10^*^ which 
currently limits the accuracy of theoretical predictions of the transition frequencies. 

Another point concerns the uncertainty of the reference data for the nuclear masses. The 
masses given in Ref. [20, 21] are the atomic masses. The mass of the bare nucleus needed for 
our calculation can be obtained by subtracting the masses of the electrons and their binding 
energies as described in Ref. ll47l . Useful tabulations of ionization energies can likewise be 
found in Ref. |!47ll48l. 

The second source of theoretical uncertainty is due to uncalculated higher-order terms 
in the theoretical expressions discussed in Sec. 12.11 The uncertainties of i-'RR and (^qed are 
estimated as follows. For the higher-order recoil terms, we use the magnitude of the last 
contribution on the right-hand side of Eq. ^ times {Za)\n[{Za)^'^]. The uncertainty in 
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Table 1. Values of the Aeo coefficient obtained with the analytical method compared to the 
numerical results for the self-energy remainder function GsE{Za) for n = 5, . . . , n = 8 
with £ = n — 1 and £ = n — 2. The numbers in parentheses ai'e standard uncertainties in 
the last figure and represent the uncertainty of the final numerical integration due to the finite 
number of lattice points (in the case of Ago), and the finite number of integration nodes for the 
numerical evaluation of multi-dimensional integrals (in the case of Gse)- 



n 


e 


2j 


K 


^60 


GsEiZ = 3) 


Gse{Z = 6) 


Gse{Z — 8) 


5 


3 


5 


3 


0.002 403 151 41(5) 


0.002 41(3) 


0.002 43(3) 


0.002 45(3) 


5 


3 


7 


-4 


0.008 087 015 45(5) 


0.008 09(3) 


0.008 11(3) 


0.008 12(3) 


5 


4 


7 


4 


0.000 814 414 71(5) 


0.000 82(3) 


0.000 82(3) 


0.000 82(3) 


5 


4 


9 


-5 


0.002 412 929 08(5) 


0.002 41(3) 


0.002 41(3) 


0.002 42(3) 


6 


4 


7 


4 


0.000 827 467 81(5) 


0.000 83(3) 


0.000 83(3) 


0.000 83(3) 


6 


4 


9 


-5 


0.002 748 250 60(5) 


0.002 75(3) 


0.002 75(3) 


0.002 75(3) 


6 


5 


9 


5 


0.000 326 676 27(5) 


0.000 32(3) 


0.000 33(3) 


0.000 33(3) 


6 


5 


11 


-6 


0.001 008 201 08(5) 


0.00101(3) 


0.00101(3) 


0.00101(3) 


7 


5 


9 


5 


0.000 325 902 82(5) 


0.000 32(3) 


0.000 33(3) 


0.000 33(3) 


7 


5 


11 


-6 


0.001 141 603 10(5) 


0.001 14(3) 


0.001 14(3) 


0.001 14(3) 


7 


6 


11 


6 


0.000 147 439 22(5) 


0.000 15(3) 


0.000 15(3) 


0.000 15(3) 


7 


6 


13 


-7 


0.000 485 185 97(5) 


0.000 48(3) 


0.000 49(3) 


0.000 49(3) 


8 


6 


11 


6 


0.000 144 496 71(5) 


0.000 13(3) 


0.000 15(3) 


0.000 14(3) 


8 


6 


13 


-7 


0.000 545 933 41(5) 


0.000 54(3) 


0.000 55(3) 


0.000 55(3) 


8 


7 


13 


7 


0.000 072 861 41(5) 


0.000 06(3) 


0.000 07(3) 


0.000 07(3) 


8 


7 


15 


-8 


0.000 258 766 38(5) 


0.000 25(3) 


0.000 26(3) 


0.000 26(3) 



^^QED is separated into two parts. For the unknown Agi term in the analytic expansion, we 
base our estimate on the analytic approach and estimate the magnitude of the ^gi term to 
be equal to Ago times (Za)'^ \n[{Za)~'^]. This provides for a uniform estimate which is 
not restricted to those states for which nonperturbative, numerical results for the self-energy 
remainder function are available. Based on a comparison of the results in states with ^ < 5, 
the magnitude of Bgo f371 is estimated to be equal to the magnitude of 4 Ago- This is used as 
the uncertainty for this second term. These estimates are analogous to those used previously 
inRef. lfT6l. 

We now illustrate, by way of example, the predictive power of QED theory for Rydberg 
state transitions in the infra-rad and near optical range, which are accessible to frequency 
comb measurements in one-electron ions. To this end, we select a two-photon transition 
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Table 2. Fundamental constants and masses used as input parameters for the evaluation of the 
theoretical expression and error estimates. In parentheses, we indicate the standard uncertainty. 
The masses m^i (A/") con'espond to the atomic mass of an atom (including the bound electrons) 
with nucleus Af. By contrast, the nuclear mass is denoted as mjv (A/") in this article (it excludes 
the mass of the bound electrons). 



Constant Value 



RooC 


3.289 841960 361(22) x lO^^ Hz 


a 


7.297 352 5376(50) x 10"^ 


ae 


1.159 652180 73(28) x 10~^ 


me 


5.485 799 0943(23) x 10-*u 




12.000 000 000(0) u 


mA("'C) 


10.016 853 2(4) u 


mA(«Li) 


6.015 122 795(16) u 


mA(^Li) 


7.016 003 4256(45) u 




1.007 825 032 07(10) u 




2.014 101 778 040(80) u 



Table 3. Theoretical predictions for two-photon transition frequencies in atomic hydrogen 
and deuteiium. The transition from the initial level |1) with quantum numbers n = 9, 
£ = 8, and j = 15/2 to the level |2} with quantum numbers n = 16, i = 10, and 
j = 19/2 in considered. For the upper state, the higher-order self-energy coefficient reads 
Aeoin = 16,e = 10, j = ^) = 1.026705(5) x 10"^. The individual contributions are 
listed in Eg. 



Term 


i/(THz) 


u{THz) 


tti 


27.749 282 6987(2) 


27.756 833 2542(2) 




0.000 000 0000 


0.000 000 0000 




0.000 000 0000 


0.000 000 0000 




0.000 000 0035 


0.000 000 0035 


Total 


27.749 282 7022(2) 


27.756 833 2577(2) 



in hydrogen, lithium and carbon and add the theoretical contributions to the transition 
frequency ([T]i given in Sec. l2.1l in order to obtain a theoretical prediction. The results are given 
in Table[3]for hydrogen and deuterium and in Table|4]for two lithium isotopes. In both cases 
we consider a two-photon transition from the state |1) with quantum numbers n — 9, £ = 8, 
j = 15/2 to a state |2) with quantum numbers n = 16, £ = 10, j = 19/2. The hydrogen- 
deuterium isotope shift with mass numbers 1 and 2 as well as for the lithium isotopes with 
mass numbers 6 and 7 is being considered. The same is done in Table |5] for another two- 
photon transition from the state |1) with quantum numbers n — 13, £ ~ 11, j = 21/2 to a 
state |2) with quantum numbers n = 17, £ = 13, j = 25/2 in carbon isotopes. If we wish to 
extract the Rydberg constant from a measurement of a two-photon transition [16..17J . then we 
have to consider ions with very well determined and known nuclear masses. In our example 
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Table 4. Theoretical predictions for two-photon transition frequencies in two isotopes of 
hydrogen-like lithium. The transition from the initial level 1 1) with quantum numbers n = 9, 
£ = 8, and j = 15/2 to the level |2) with quantum numbers n = 16, £ = 10, and j = 19/2 
in considered. We recall that for the upper state, the higher-order self-energy coefficient reads 
Aeoin = 16,£ = 10, j = ^) = 1.026705(5) x 10~^. Again, the individual contributions 
are listed in Eq. (9). 



Term 


6Li2+ i/(THz) 


^Li2+ !/(THz) 


I^D 


249.857 322 9816(17) 


249.860 575 0963(17) 


i^G 


0.000 000 0000 


0.000 000 0000 




0.000 000 0000 


0.000 000 0000 


'^QED 


0.000 000 2829 


0.000 000 2829 


Total 


249.857 323 2645(17) 


249.860 575 3793(17) 



Table 5. Transition frequencies for the transition from level |1) with quantum numbers n = 
13,£ = ll,andj = 21/2to level |2) with quantum numbers n = 17, £ = 13, and j = 25/2 
in two isotopes of hydrogen-like carbon. The self-energy remainder for the upper state is 
estimated based on a coefficient of A6o(n = 17, £ = 13, j = ^) = 3.76900(5) X 10"^. 



Term 


I2c5+ KTHz) 


I0c5+ i/(THz) 


I'D 


290.976 045 2425(19) 


290.973 410 2736(20) 


i'BG 


0.000 000 0000 


0.000 000 0000 


i^R 


0.000 000 0001 


0.000 000 0001 


''QED 


0.000 000 6216 


0.000 000 6216 


Total 


290.976 045 8641(19) 


290.973 410 8953(20) 



cases considered above in Table [5] the relative accuracy of the electron to nucleus mass ratio 
is as follows: 

r(iH) - ' r(2H) ' ^.i x 10 . (15) 

By contrast, the relative accuracy of the nuclear masses for the example cases in Tables |4] 
and|5]is as follows, 

Currently, none of the given uncertainties in the mass ratios limit the final accuracy of the 
theoretical predictions of Rydberg state transitions (the theoretical accuracy is limited by the 
current value of the Rydberg constant, on a level of about 7 x lO^^^). However, in view 
of a conceivable improvement of the accuracy of the Rydberg constant in the near future, 
the electron to nucleus mass ratios of and as well as the nuclear masses of ^Li, ^Li 
and may soon become a limiting factor. Specifically, for frequency measurements better 



Nuclear masses and spectroscopy of Rydberg states 



11 



than 6.0 x 10^^'* for ^Li, the nuclear mass accuracy will become limiting, and for the other 
nuclei, this effect becomes relevant at a relative accuracy of about 10^^^ for the frequency 
measurements. Conversely, frequency measurements of better accuracy allow for a better 
mass determination using spectroscopy. The latter aspect will be analyzed in more detail in 
the following. 

3. Nuclear Mass Determination 

3.1. General Paradigms 

We intend to explore the applicability of three methods for the determination of nuclear 
masses via high-precision spectroscopy of Rydberg transitions in hydrogen-like ions. As 
already outlined in Sec. [T] one of these methods is essentially based on a very precise 
measurement of the isotope shift of a specific Rydberg transition for two isotopes, one of 
which has a very well known reference mass. We denote the nucleus with the accurately 
known reference mass as Afn and its mass as tojv(7V/;,), whereas the other nucleus, whose 
mass mN{MM) is to be determined, is denoted as Mm- One measures two transition 
frequencies, one in the reference system iv^j^^ and one in the isotope whose mass is to 
be determined (denoted v^l^2)- The system of equations composed of the two measured 
transition frequencies v^^2 ^"d Vii^2 ^^'^ QED theory can then be solved for the Rydberg 
constant and for the mass of the unknown isotope. This method is known as method I 'm the 
current paper 

For method II and method III, we assume that various efforts of deducing an improved 
value for the Rydberg constant | T6l [TTl |26J using a Rydberg transition in a hydrogen-like 
ion with a very well known nuclear mass are successful. A second measured transition 
frequency !^f^2 ^ system with an inaccurately known nuclear mass mM{MM) can directly 
be compared to its theoretical value. The nuclear mass mN{J\fM) or the electron to nucleus 
mass ratio r{MM) can then be determined. 

For the above methods, the following order-of-magnitude estimates are relevant: We lose 
about four decimals when converting the frequency measurement to a measurement of the 
nuclear mass. Based on a comparison to high-precision spectroscopy of lower-lying levels 
of atomic hydrogen f^^l, we assume that a reasonable target accuracy for a Rydberg state 
transition lies in the range of 10^ A sufficient accuracy of a nuclear reference mass for our 
purposes therefore corresponds to the level of 10"^*^. Two Rydberg transition frequencies of 
relative accuracy 10^^"*, one of which in a reference system, can be solved for the Rydberg 
constant (of accuracy 10^^"*) and for the unknown nuclear mass (to be determined with an 
accuracy of 10^^*^, in the context of method I). Alternatively, if the Rydberg constant is 
independently known with an accuracy of 10^^'' via a measurement in a reference system, 
then any other nuclear mass can be determined with an accuracy of 10^^° by a measurement 
of a different transition in a different hydrogen-like ion (the latter consideration is relevant to 
method 11). 

Let us now cast these order-of-magnitudes estimates into formulas. We first have to 
express the dependence of a transition frequency i^i<^2 on the electron to nucleus mass ratio 
r(7V) = me/mN{Af), where tojv(7V) is the mass of the bare nucleus (without the bound 
electrons). As can be seen from the theoretical expressions in Sec. 12.11 this dependence has 
a complicated functional form. E.g., in the Dirac value Eq. (|2]l, the first term is directly 
proportional to the reduced mass fir = 1/(1 + ^{M")), while the second is proportional to 
r{Af)iX^- However, the dominant and leading dependence is simply given by the approximate 
proportionality of the transition frequency to the reduced mass of the system. We can thus 
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define a scaled frequency fi^2^ which is related to 1^14^2 by 

1 



J^l<^2 — RooC 



■h 



(18) 



1 + r{J\f) 

The scaled frequency fi^2 given by theory still carries a residual dependence on the mass 
ratio. For stable/long living nuclei equal in mass or heavier than lithium, known electron to 
nucleus mass ratios provide enough accuracy so that the residual dependence of fi^2 does 
not contribute to the uncertainty on a level required for the nuclear mass determination. 

Formulas for method I: We need one isotope with a well-known nuclear mass which we 
will use as the reference system. Its transition frequency, mass ratio and theoretical value are 
labeled Vij^2^ fV^R)^ and //^2' respectively. Based on this reference, we want to measure 
the mass of the nucleus of another isotope. Its transition frequency, mass ratio and theoretical 
value will be labeled Vii^2^ ^{Mm), and fil^2- measuring the transition frequencies Vij^2 
and vf^l^2 in both isotopes, we get the system of equations 

RooC , V . . . (19) 



1 



M 

a«2 



1 



r 



[Um] 



M 
«2 



(20) 



We cancel the Rydberg constant and obtain 



,M 
a-H-2 



M 
02 



mNiMn) f 1^2 mNjAfAi) + 



rUNiJ^R) +me 

Solving for the nuclear mass mN{NM), we obtain 



(21) 



(22) 



X TOg 



\'^2 



M 
i-i-2 



,M 



R 

<-»-2 



{rfiNiMR) + rrie) - mN{NR) 



This allows us to determine the nuclear mass of one isotope mNiMM) from a measurement 
of a transition frequency J^f^2 this isotope, and a reference transition frequency v^^2 an 
isotope with nuclear mass mpf{AfR). 

Formulas for method 11: We have already mentioned a joint theoretical and experimental 
project with the National Institute of Standards and Technology (NIST), whose aim is to 
deduce an improved value for the Rydberg constant from transitions in Rydberg states lfT6l 
\T7]. Furthermore, we also mention a project at the National Physics Laboratory (NPL) in the 
United Kingdom where the 2S-8D transition in hydrogen is intended to be used in order to 
improve the accuracy of the Rydberg constant ll26l . Let us assume that one of these efforts is 
met with success and the uncertainty of the Rydberg constant can be reduced significantly (by 
at least an order of magnitude as compared to the 2006 COD ATA value iTT\ ). 

This would allow to open up another possible way of measuring masses which is 
especially interesting for isotope systems where no mass is known well enough to serve as a 
reference mass. Namely, provided the Rydberg constant can be determined to good accuracy 
in a different ionic system, the nuclear mass mM{MM) can be obtained by solving Eq. ( fTSl l 
yielding 

' fM n \ ^1 

J1^2^o 



TUN {Ml 



M) 



,M 



1 



(23) 
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Numerical loss is incurred because 

■^^^L -1=KA/'m), (24) 

where r{MM) is rather small 10~^ . . . 10~^ in typical cases), whereas the two terms on 
the left hand side are of order unity. 

This equation also allows for a determination of the electron to nucleus mass ratios which 
we will denote as method III in the following. In this case, we just use 

t{Nm)=^-^^^P^-1. (25) 

in order to determine r{MM), again with a loss in numerical significance of above four 
decimals. 

Table 6. We list the sources for the relative uncertainties for the determination of 
(5m :\r (A/jvf ) /mjv (A^M ) We explore method I for the determination of a nuclear mass (see 
text). The reference nucleus is A^^ = -'^•^C, and the nuclear mass of A^a/ = ^"C is to be 
measured. The transition is from state |1) <r-> |2), where |1) has quantum numbers ii = 1,3, 
£ = 11 and J = 21/2. whereas |2) has quantum numbers n = 17, £ = 13 and j = 2-5/2. The 
contributions to the relative uncertainty (5m jv (A^M ) /m]v (A/m ) above the horizontal line are 
evaluated in terms of the 2006 CODATA recommended values of the fundamental constants 
and do not influence the determination of mff(J^M) on the level of one part in 10"^°. The 
contributions below the horizontal line are due to the assumed experimental (spectroscopic) 
accuracy of the transitions. 



Source 



Stun {.Mm) 
TUN (Mm) 



1 dmMjJ^M) 

mN{MM) dmN(MR) 

1 dniNiAfM) 
mNWrn) dr{MM) 

1 dmN{MM) 
mN{MM) dr{MR) 

1 druNiJ^M) 
ruNiNM) da 

1 amjv(A^M) 

'"■^(A/'m) Sle 



SniMiNa) 9.6 X 10~" 

6r(NM) 1.8 X 10-" 

5r(MR) 1.6 X 10-15 

(5a 1.0 X 10-15 

Sae 8.8 X 10-20 



1 druNiMM) ; „„ 
T^^^ -Srrie 6.9x10 

S[f^^2/fS^2] 1-6 X 10-" 

^JI^EpA 5u^^, 1.8 X 10* f ^ 



.2 



l<->-2 



3.2. Potential of the Method 

We now illustrate the two discussed methods for the determination of nuclear masses from 
Rydberg transition spectroscopy in hydrogen-Uke ions. To this end, we keep the attained 
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Table 7. We explore the application of method II (see text) for the determination of the mass 
of the A/a/ = ^Li nucleus. The transition is |1) |2) where |1) is the state with quantum 
numbers n = 9, £ = 8, and j = 15/2, and |2) has quantum numbers n = 16, £ = 10 
and j = 19/2. Contributions because of theoretical input data to the relative uncertainty 
Sm^ {Afj^.j) /mpf {J\f]\{) of the nuclear mass are given above the horizontal line, whereas 
contributions due to the assumed spectroscopic measurements are given below the horizontal 
Une. 



Source 



mN(XM) 



1 diUNiAfM) 







1 


druMiAfM) 




da 


1 


dmN{NM) 




dae 


1 


9mjv(A/'M) 



5r{MM) 1-5 X 10-15 

5a 4.4 X 10-" 

(5ae 3.6 X 10-15 

5m, 4.2 X 10-1° 



&flL2 1-1 X 10- 

1 dmM{MM) rM 1 Q -,n4 ( 



5uf^2 1.3 X 10^ 



•f+2 \ ''l<->2 



1«2 



1 drriNiJ^M) irt^mif^^^ 
■ oHooC 1.3 X 10 



experimental accuracy as a variable and show to which accuracy nuclear masses can be 
determined for a given assumed spectroscopic accuracy. We concentrate on transitions in 
isotopes of atomic hydrogen, lithium and carbon that lie in the infra-red frequency range, 
which is an ideal appUcation range for frequency combs. 

For this illustration, we use the atomic mass unit which is also the unit generally used 
in measurements in Penning traps. There the masses are determined by a comparison of 
cyclotron frequencies of different ions. Most atomic masses and especially the mass of the 
electron which is rather important for our methods are in fact known more precisely in the 
atomic mass system. Even though there are at the moment numerous efforts to use the atomic 
mass unit to define the SI unit kilogram, so far the conversion factor between the two units 
still carries a rather large relative uncertainty of 5.0 x 10~^. In the determination of mass 
ratios (close to unity) or mass differences, the uncertainty of the conversion factors become 
irrelevant, and these are the quantities that one needs for applications to fundamental physics. 
We can thus safely work in atomic mass units. 

We first discuss an example application for method I, which as we recall implies the 
determination of Rydberg constant and nuclear mass from an accurate measurement of two 
transitions in hydrogen-like ions of different isotopes. The reference system is taken as the 
one-electron ion of ^^C. We recall that the ^^C atom serves as the definition of the atomic 
mass unit (u) which is defined by the relation TOyi(^^C) = 12 u. The nuclear mass tojv(^^C) 
is obtained from to^(^^C) by the subtraction of the electron rest masses and the binding 
energies. The latter carries a theoretical uncertainty which has been discussed in [47J; yet 
mAr(^^C) still represents a very accurate nuclear mass standard. Even taking into account 
this uncertainty, the final relative uncertainty of the nuclear mass is 1.2 x 10~^^ which is 
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Table 8. We explore the application of method III (see text) for the determination of the 
electron to nucleus mass ratio for A/m = ^H. The transition is |1) |2) where |1) is the state 
with quantum numbers n = 9, £ = 8, and j = 15/2, and |2) has quantum numbers n = 16, 
(, = 10 and j = 19/2. The contributions to the relative uncertainty 5r{MM)/r{MM) 
of the electron to deuteron mass ratio due to the 2006 CODATA values for the fundamental 
constants and masses required for the evaluation of the theoretical expressions are given above 
the horizontal Une. Contributions due to the assumed spectroscopic measurements are given 
below the horizontal Hne. 



Source 



r{NM) 



1 








1 


dr(MM) 


r{MM) 


da 


1 


dr{NM) 



5r(NM) 9.1 X 10-" 

5a 1.4 X 10-12 

5ae 1.2 X 10-18 



r(MM) dae 
5f^^2 6.8 X 10-33 

1 dri^Ml M 3 7 X 103 r 



1 dr{NM) 
r(NM) dRooC 



5RooC 3.7 X 10^ 



V RooC ) 



Table 9. We explore the application of method III (see text) for the determination of the 
electron to nucleus mass ratio for A/m = ^H. The transition is |1) |2) where |1) is the state 
with quantum numbers n = 9, £ = 8, and j = 15/2, and |2) has quantum numbers n = 16, 
^ = 10 and j = 19/2. Again, contributions to the relative uncertainty Sr{MM)/r{MM) 
of the mass ratio are separated into those caused by theoretical input data which are given 
above the horizontal line, and contributions due to the assumed spectroscopic measurements 
are given below the horizontal line. 



Source 



&r{NM) 
r{MM) 



1 drjAfM) 
r(AfM) dr{I^M) 

1 dr(NM) 
''{J^m) da 

1 driNM) 



5r{AfM) 9.3 X 10-" 

5a 7.0 X 10-" 

5ae 5.8 X 10-" 



r(A/'M) dae 
5f^^2 1-3 X 10-32 



1 dr(MM) ^^.M isxlO^^f^^ 

/ 5R^c \ 
V RooC ) 



1 dr{NM) 
r{NM) dRooC 



5RooC 1.8 X 103 



Still two orders of magnitude more accurate than any measured nuclear mass. The isotope 
whose nuclear mass we want to determine is assumed to be the one-electron ion of ^°C. In 
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both isotopes, the same transition is employed for the mass determination. The currently best 
value for the atomic mass of is given in Table|2]and has a relative uncertainty of 4 x 10^^. 
Although this mass determination to 4 parts in 10^ is very accurate on an absolute scale, we 
should remember that there appears to be room for improvement. For comparison, one of the 
most precise mass measurements has recently been conducted for ^^O with a relative precision 
of 1.1 X 10-" (seeRef. ISOl ). 

In terms of our notation for the "master equation" Eq. (l22l i for method I, we have 
Mr = ^^C and Mm = ^"C. We specify |1) ^ |2) as the same transition whose frequency 
(based on the current values of the fundamental constants |ir|) has akeady been evaluated 
in Table |5] Above the horizontal line of Table |6] we list the contributions to the relative 
uncertainty 5mAr (A/'m ) at (A/jvf ) of the measured mass of the Mm = due to the 
reference nucleus mass, due to the electron-to-nuclear mass ratios r{AfM) and ^(-Vr), due 
to the 2006 COD ATA value of the fine-structure constant, due to the experimental uncertainty 
B6l of the electron magnetic moment anomaly a,, which enters the theoretical expression in 
Eq. and due to the higher-order QED terms which enter the theoretical expressions for the 
scaled transition frequencies ff!^2 ™d //^2- chosen transition has j = t — \ for both 
states involved. In general, states with j = t — \ allow for QED predictions with slightly 
smaller uncertainties than those with J = ^ + ^. This is because states with j = ^ — ^ in 
general have numerically smaller coefficients in the semi-analytic representation of radiative, 
radiative-recoil and recoil corrections as compared to states with j = ^ + |- This trend is 
expected to continue for the uncalculated higher-order terms. Our procedure for estimating the 
higher-order terms, which expresses the estimates in terms of multiples of known coefficients, 
therefore leads to slightly smaller theoretical uncertainties for the mentioned transitions as 
compared to a case with j = £ + i. 

If we assume that the experimental accuracy for the spectroscopic measurements of 
5vf^2 and 5v^l^2 is of the order of lO"^"*, then the by far dominant contribution to the 
uncertainty of such a mass determination is caused by the uncertainty in the frequency 
measurements. Their contribution to 5mM{.N'M) /"mNi-I^M) is roughly four orders of 
magnitude larger than Svij^2l^u^2 ™d 5vfl^2l'^iU2^ as given below the horizontal line 
in Table |6] Finally, to give a numerical example, we conclude that assuming a relative 
uncertainty in the frequency measurements of the order of 1.4 x 10"^^ could be reached, 
the nuclear mass of ^'^C could be determined with a relative uncertainty of 3.6 x lO"* on a 
level comparable to the 2006 CODATA value. We hereby add all contributing uncertainties 
quadratically. For completeness, we note that even if QED theory carried no uncertainty, 
and even if our spectroscopic accuracy were perfect, the maximum accuracy reachable with 
method I would still be limited by the uncertainty of the reference mass. Therefore having a 
reference mass with a high accuracy such as ^^C is recommendable in order to minimize one 
particularly problematic source of uncertainty. 

We now discuss the application of method II, which implies a measurement of a single, 
specific transition in a hydrogen-like ion, with the Rydberg constant having been determined 
in a different system. The "master equation" for this method is Eq. (l23T l. Here, the example 
we study is Mm = ^Li, with the transition given in Table |4] The current best value |20| for 
the mass of ^Li has a relative uncertainty of 6.4 x IQ-^". The difference, however, to the 
previously accepted value for the mass of ^Li 121] is 14 standard deviations or 1.1 x IQ-^u. 
Another measurement with an independent method with a comparable accuracy to f2n\ would 
therefore be of great interest. In Table|7] we list the different sources of theoretical uncertainty 
for the mass determination. As evident from the entry in the fifth row of Table|7] the currently 
accepted value of the electron mass also is significant for the mass determination of the lithium 
isotope, though not the primary limiting factor. It is generally assumed that an improved value 
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of the electron mass could be within reach of an improved measurement of a bound-electron 
g factor in a low-Z hydrogen-Uke ion ||5T] l52l |53l |5?1 1551 . 

If it would be possible to reach an accuracy of the frequency measurements in the 
range of 2 x 10"^'* for both the Rydberg determination frequency as well as the Lithium ion 
frequency (in the setting discussed above), then the nuclear mass of ^Li could be determined 
with a relative uncertainty of 5.6 x lO^^*'. Thereby, a similar accuracy as in the measurement 
of II20I could be reached. 

In order to examine whether the required precision in the measurement can be reached 
for the discussed transitions, the ratio of the transition energy E to the decay width of the line 
r has to be considered. In lfT6l . based on a calculation for the decay width for Rydberg states 
carried out in ll56l . this Q factor has been evaluated for transitions of near circular Rydberg 
states from n to n — 1 as 

En - En-i 3n^ 

For the transitions considered here a complete calculation of the decay width following the 
description in [57 1 leads to the results for the Q factor 

Q(Li)[n = 12,f 10 n = 9,^ = 8] = 8.0 x 10^ (27) 
Q(H)[n = 12,^= 10^n==9,^ = 8] = 7.2 x 10^ , (28) 

Q(C)[n = 17,^ = 13 ^ n = 13,£= 11] = 2.8 x 10^ (29) 

The comparison of the Q factors for hydrogen and lithium show its dependence on the nuclear 
charge number Z. Measurements in the past have been able to determine the energy of a 
transition within lO^"* of the width of the line |8 58 |. It would be required to enhance this to 
about 10^^ to reach the required accuracy for a mass determination in Li and to 3 x 10^ for a 
mass determination in C. The attempt at NIST fT6',T7'| to measure transitions among Rydberg 
states up to an accuracy of 10^^^ might requires one to develop techniques which enhance 
the resolution of the lines to the indicated values. 

Due to the higher Q values for hydrogen, the uncertainty with which the transition 
frequencies can be determine would be about one order of magnitude smaller Moreover, 
the electron to proton or deuteron mass ratios are much larger compared to mass ratios for 
nuclei with higher Z considered so far Because the numerical loss is directly related to the 
mass ratio [see Eq. (l24bl. another source of uncertainty is reduced. 

This leads us to the discussion of the application of method III. Here, we investigate the 
two cases Mm — and Mm = using the transitions in Tab. [3] with the master equation 
(I25] ) of method III. The theoretical sources of uncertainty are given in Table [8] Based on 
Refs. |I8]|58|, we assume that the line can be split to at least one part in 10"^. Together with 
the Q factor this leads to a conservative estimate of the relative accuracy of the frequency 
measurement of 1.4 x 10"^'^, which we will also use for the assumed relative uncertainty of 
the Rydberg constant. These uncertainties would allow to determine the electron to deuteron 
mass ratio with a relative uncertainty of 7.3 x 10^^°; this is comparable to the 2006 CODATA 
mass ratio which has a relative uncertainty of 4.2 x 10^^° ifTTI . Every improvement in the 
resolution of the line would allow to increase the accuracy with which the mass ratio can be 
determined. For A/a/ = ^H, the theoretical sources of uncertainty are listed in Table |9l We 
use the same conservative estimate as for regarding the uncertainty of the measurement of 
the transition frequency and the Rydberg constant. Even so, the electron to proton mass ratio 
could still be determined with a relative uncertainty of 3.6 x 10^^°; this would be slightly 
better than the current best value which has a relative uncertainty of 4.3 x 10" ^'^ li2U . 
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4. CONCLUSIONS 

The primary goal of the current work has been to show that high-precision spectroscopy of 
Rydberg states in hydrogen-like ions has a far greater potential for advances in fundamental 
physics than just the determination of the Rydberg constant and can be used for the 
determination of nuclear masses as well. This finding enhances the likely impact of a 
successful implementation of high-precision spectroscopy of highly excited Rydberg states 
of hydrogen-like ions with a low and medium nuclear charge number In general, we can say 
that the basic idea of the current paper is to use a lightly bound electron (in a Rydberg state) to 
probe the nucleus (a "Rydberg electron trap" for the nucleus). The interaction of the nucleus 
to the loosely bound electron could be measured with high precision. Because this interaction 
is well understood, it can be used to determine the mass of the nucleus. Hyperfine effects 
caused by the interaction of the nuclear spin with the total angular momentum of the electron 
are briefly considered here in Eq. If necessessary, they can be calculated to relative order 
a{Za)^ suing the formalism recently outlined in Ref. |41 1. 

The theoretical predictions for the transition frequencies in Rydberg states are more 
precise than for lower-lying transitions, because a number of problematic effects like the 
nuclear-size correction are effectively suppressed. It is thus possible to devise transitions 
in the infra-red (see Tables [3]— for probing nuclear masses; these transitions combine a 
favourable range of frequencies for the application of optical frequency combs with small 
uncertainties in the theoretical QED predictions. Summarizing the theoretical calculations, we 
can say that with an entirely realistic relative uncertainty of 1.4 x 10^^^ for the measurement 
of the transitions in deuterium and hydrogen, the electron to proton mass ratio could be 
determined with a relative uncertainty of 3.6 x 10^^". This would mean a small improvement 
over the present value with a relative uncertainty of 4.3 x 10^^''. For the electron to deuteron 
mass ratio, a relative uncertainty comparable to the 2006 CODATA value could be achieved. 
If it were possible to improve the accuracy further and to reach a relative uncertainty of 
2 x 10^^"* for the two-photon transition in hydrogen-like ^Li, we could improve lithium 
nuclear masses by more than one order of magnitude as compared to their values in the 
AME2003 [21 1 adjustment and of similar accuracy as 1 20J . This would potentially confirm 
the large difference between these two values for the ^Li mass from Refs. Il20l 1211 by an 
independent method. As evident from the entry in the fifth row of Table |7] the accuracy 
could be further improved with a more precise value for the electron mass whose accurate 
determination is currently being pursued, e.g., in precise measurements of the bound-electron 
g factors. The mass determination of using the isotope shift would allow us to reach a 
relative uncertainty comparable to the current relative uncertainty of its mass of 4.0 x 10^^ 
II2TI if the transition frequency could be measured with a relative uncertainty of 1.4 x 10^^^. 

As evident from the last two rows of Tables |6]|7] |8]and|9] one loses about four decimals 
in the determination of a nuclear mass from transition frequencies, due to inevitable numerical 
cancellations. For method II, this is explained in the text following Eq. (l24l l. Nevertheless, 
the approach leads to very promising results. With the progress made in recent years in 
reducing the relative uncertainties in the frequency measurement using frequency combs (see 
||59l and references therein), there is justified hope for even more precise results for transition 
frequencies in the future. In turn, they would allow to further increase the accuracy of such a 
mass determination. In view of the present efforts of measuring nuclear masses with increased 
accuracy, it certainly would be helpful to have an independent method to check the nuclear 
masses determined by comparing cyclotron frequencies in Penning traps (60]. 

Concerning the broader impact of the current proposal, we remember that ^Li is one 
of the primordial elements, and a more precise mass measurement may help to improve our 
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understanding of the processes in the Big Bang Nucleosynthesis, as well as the reason for 
the low binding energy of the discussed nucleus. Moreover, ^Li also plays a role in many 
nuclear reaction processes especially in the fusion cycle of the sun where an improved mass 
measurement can help to find bounds for the mass of the solar neutrino |19| and also for 
the determination of other masses through nuclear reactions ||6T1 . The electron to proton 
and electron to deuteron mass ratios are important for many applications in spectroscopy 
where these two isotopes are studied very intensively. For example, a more accurate electron 
to proton mass ratio might be very interesting for the planned comparison of transition 
frequencies in hydrogen and anti-hydrogen ifTSl l62l in order to test CPT invariance. Both 
mass ratios also play a role in the analysis of the hydrogen-deuterium isotope shift which 
is used to deduce the difference of the mean square charge radii of the proton and the 
deuteron [3]. Another field where precise nuclear masses are of great interest is the study 
of the unitarity of the CKM matrix 1631 l22l . where super-allowed beta decays are being 
investigated to determine the CKM matrix element Vud- The nucleus is of particular 
interest in these studies |23|. A further increase in accuracy for the nuclear mass tojv(^°C) 
can improve the bounds on CKM non-unitarity and therefore provide a check for the standard 
model. 

In order to provide reliable theoretical predictions for such joint experimental-theoretical 
efforts, we have performed analytic as well as numerical calculations of the self-energy 
remainder function for states with n = 5, . . . , 8 with £ = n — 1 and £ = n — 2, and we 
have also obtained values for the Aqq coefficient for these states. The theoretical predictions 
for the notoriously problematic self-energy remainder functions show an excellent mutual 
agreement. This leaves little room for any conceivable changes in the theoretical predictions 
due to calculational errors whose existence in complex QED calculations otherwise cannot 
be ruled out without extensive cross checks. The field appears to be open for experimental 
studies. 

Acknowledgments 

Insightful discussions with P. J. Mohr and J. N. Tan are gratefully acknowledged. This 
research has been supported by the National Science Foundation (Grant PHY- 8555454) and 
by a Precision Measurement Grant from the National Institute of Standards and Technology. 

References 

[1] Niering M, Holzwarth R, Reichert J, Pokasov P, Th Udem, Weitz M, Hansch T W, Lemonde P, Santarelli G, 
Abgrall M, Laurent R Salomon C and Clairon A 2000 Phys. Rev. Lett. 84 5496-5499 

[2] Weitz M, Huber A, Schmidt-Kaler F, Leibfried D, Vassen W, Zimmermann C, Pachucki K, Hansch T W, Julien 
L and Biraben F 1995 Phys. Rev. A 52 2664-2681 

[3] Huber A, Th Udem, Gross B, Reichert J, Kourogi M, Pachucki K, Weitz M and Hansch T W 1998 Phys. Rev. 
Leff. 80 468^71 

[4] de Beauvoir B, Nez F, Julien L, Cagnac B, Biraben F, Touahii D, Hilico L, Acef O, Clairon A and Zondy J J 

1997 Phy.s: Rev Lett. 78 440^43 
[5] Schwob C, Jozefowski L, de Beauvoir B, Hilico L, Nez F, Julien L, Biraben F, Acef O, Clairon A and Zondy 

J J 1999 Phys. Rev. Lett. 82 4960^963 
[6] de Beauvoir B, Nez F, Julien L, Cagnac B, Biraben F, Touahii D, Hilico L, Acef O, Clairon A and Zondy J J 

1997 Phys. Rev Lett. 78 440^43 
[7] Berkeland D J, Hinds E A and Boshier M G 1995 Phys. Rev Lett. 75 2470-2473 
[8] Hagley E W and Pipkin F M 1994 Phys. Rev Lett. 72 1 172 
[9] Lundeen S R and Pipkin F Ml 986 Metrologia 22 9 
[10] Newton G, Andrews D A and Unsworth P J 1979 Phil. Trans. R. Soc. A 290 373 
[1 1] Mohr P J, Taylor B N and Newell D B 2008 Rev. Mod. Phys. 80 633 



Nuclear masses and spectroscopy ofRydberg states 



20 



[12] Pachucki K 1996 Phys. Rev. A 53 2092 
[13] Pachucki K 1999 Phys. Rev. A 60 3593 

[14] Friar J L, Martorell J and Sprung D W L 1997 Phys. Rev. A 56 4579^586 

[15] Jentschura U D, Kotochigova S, Le Bigot E O, Mohr P J and Taylor B N 2005 Phys. Rev. Lett. 95 163003 
[16] Jentschura U D, Mohr P J, Tan J N and Wundt B J 2008 Phys. Rev. Lett. 100 160404 
[17] Jentschura U D. Mohr P J, Tan J N and Wundt B J 2009 Can. .J. Phy.s. 87 757-762 

[18] Gabrielse G, Larochelle P, Le Sage D, Levitt B, Kolthammer W S, McConnell R, Richerme P, Wrubel J, Speck 
A, George M C, Grzonka D, Oelert W, Sefzick T, Zhang Z, Carew A, Comeau D, Hessels E A, Storry C H, 
Weel M and Walz J 2008 Phys. Rev. Lett. 100 113001 

[19] Reines F 1967 Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences 301 

159-170 

[20] Nagy S, Fritioff T, Suhonen M, Schuch R, Blaum K, Bjorkhage M and Bergstrom I 2006 Phys. Rev. Lett. 96 
163004 

[21] Audi G, Wapstra A H and Thibault C 2003 Nuclear Physics A 729 337 - 676 

[22] Blaum K 2006 Physics Reports 425 1 - 78 

[23] Hardy J C and Towner I S 2006 Nucl. Phys. News 16N4 1 1-17 

[24] Koelemeij J C J, Roth B, Wicht A, Emsting I and Schiller S 2007 Phys. Rev. Lett. 98 173002 

[25] Karr J P, Bielsa F, Doufflet A, Pedregosa Gutierrez J, Korobov V 1 and Hffico L 2008 Phys. Rev. A 77 063410 

[26] Howers J, Baird P, Klein H, Langham C, Margolis H and Walton B 2008 Conference on Precision 

Electromagnetic Measurements Digest, 2008. CPEM2008. 40-41 
[27] Kleppner D (private communicaton) 
[28] Vries J C D Ph.D. thesis, MIT, 2001 

[29] Eides M I, Grotch H and Shelyuto V A 2001 Phys. Rep. 342 63-261 

[30] Sapirstein J and Yennie D R 1990 Theory of hydrogenic bound states Quantum Electrodynamics {Advanced 
Series on Directions in High Energy Physics vol 7) ed Kinoshita T (Singapore: World Scientific) pp 560- 

672 

[31] Barker W A and Glover F N 1955 Phys. Rev. 99 317 
[32] Erickson G W 1977 /. P/i>'.v. Chem. Ref Data 6 831 
[33] Pachucki K and Grotch H 1995 Phys. Rev. A 51 1854 

[34] Golosov E, Elkhovskii A S, Milshtein A 1 and Khriplovich I B 1995 Zh Eksp. Teon Fiz. 107 393 [JETP 80 (2), 
208 (1995)] 

[35] Jentschura U and Pachucki K 1996 Phys. Rev. A 54 1853-1861 

[36] Wichmann E H and Kroll N M 1956 P/i>'.v. Rev 101 843 

[37] Jentschura U D 2006 Phys. Rev A 74 062517 

[38] Wundt B and Jentschura U D 2008 Phys. Lett. B 659 571-575 

[39] Jentschura U D and Yerokhin V A 2006 Phys. Rev. A 73 062503 

[40] Brodsky S J and Parsons R G 1967 Phys. Rev. 163 134-146 

[41] Jentschura U D and Yerokhin V A 2010 Phys. Rev. A 81 012503 

[42] Caswell W E and Lepage G P 1986 Phys. Lett. B 167 437 

[43] Jentschura U D, Czamecki A and Pachucki K 2005 Phys. Rev. A 72 062102 

[44] Salomonson S and Oster P 1989 Phys. Rev. A 40 5559 

[45] Jentschura U D, Mohr P J and Tan J N 2009 J. Phys. B in press 

[46] Hannecke D, Fogwell S and Gabrielse G 2008 Phys. Rev. Utt. 100 120801 

[47] Mohr P J and Taylor B N 2005 Rev. Mod. Phys. 77 1 

[48] Lide D R 2007 CRC Handbook of Chemistry and Physics, 88th Edition ( Crc Handbook of Chemistry and 

Phy.sics) (Boca Raton, FL: CRC) 
[49] Fischer M, Kolachevsky N, Zimmermann M, Holzwarth R, Th Udem, Hansch T W, Abgrall M, Griinert J, 

Maksimovic 1, Bize S, Marion H, Pereira Dos Santos F, Lemonde P, SantareUi G, Laurent P, Clairon A, 

Salomon C, Haas M, Jentschura U D and Keitel C H 2004 Phys. Rev. Utt. 92 230802 
[50] Dyck R S V, Jr, Pinegar D B, Liew S V and Zafonte S L 2006 International Journal of Mass Spectrometry 251 

231 -242 

[51] Verdii J, Djekic S, Stahl S, Valenzuela T, Vogel M, Werth G, Beier T, Kluge H J and Quint W 2004 Phys. Rev. 
Lett. 92 093002 

[52] Vogel M, Alonso J, Djekic S, Kluge H J, Quint W, Stahl S, Verdu J and Werth G 2005 Nuclear Instruments and 
Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 235 7-16 

[53] Werth G, Alonso J, Beier T, Blaum K, Djekic S, Hffher H, Hermanspahn N, Quint W, Stahl S, Verd J, Valenzuela 
T and Vogel M 2006 International Journal of Mass Spectrometry 251 152 - 158 

[54] Quint W, Nikoobakht B and Jentschura U D 2008 Pis 'ma v. Zh £ksp. Teor. Fiz. 87 36^0 [lETP Lett. 87, 30 
(2008)] 

[55] Quint W, Moskovkhin D L, Shabaev V M and Vogel M 2008 Phys. Rev. A 78 032517 

[56] Jentschura U D, Le Bigot E O, Evers J, Mohr P J and Keitel C H 2005 J. Phys. B 38 S97-S105 



Nuclear masses and spectroscopy of Rydberg states 



21 



[57] Bethe H A and Salpeter E E 1957 Quantum Mechanics of One- and Two-Electron Atoms (Berlin: Springer) 

[58] Lundeen S R and Pipkin F M 1981 Phys. Rev. Lett. 46 232 

[59] Haensch T W 2006 Rev. Mod. Phys. 78 1297 

[60] Brown L S and Gabrielse G 1986 Rev. Mod. Phys. 58 233-31 1 

[61] Ball G C, Davies W G, Forster J S and Hardy J C 1972 Phys. Rev Lett 28 1069-1071 

[62] Bluhm R, Kostelecky V A and Russell N 1999 Phys. Rev Lett 82 2254-22 57 

[63] Abele H et al. 2004 Eur. Phys. J. C33 1-8 (Preprint he-p-ph / Q 312150) 



